Laplacian growth with separately controlled noise and anisotropy 
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Conformal mapping models are used to study competition of noise and anisotropy in Laplacian 
growth. For that, a new family of models is introduced with the noise level and directional anisotropy 
controlled independently. Fractalization is observed in both anisotropic growth and the growth with 
varying noise. Fractal dimension is determined from cluster size scaling with its area. For isotropic 
growth d — 1.7, both at high and low noise. For anisotropic growth with reduced noise the dimension 
can be as low as d — 1.5 and apparently is not universal. Also, we study fluctuations of particle 
areas and observe, in agreement with previous studies, that exceptionally large particles may appear 
during the growth, leading to pathologically irregular clusters. This difficulty is circumvented by 
using an acceptance window for particle areas. 

PACS numbers: 61.43.Hv, 47.54.+r, 81.10.Aj 



I. INTRODUCTION 

In a large class of pattern-forming systems the growth 
is controlled by a Laplacian field. In diffusion limited 
aggregation (DLA) this field is the probability density 
of aggregating particles fl],[§. ln viscous fingering it is 
pressure |jj , and in crystal growth it can be either a dif- 
fusive or a thermal field [Q. After the DLA model was 
introduced by Witten and Sander [||, it became stan- 
dard to simulate the Laplacian field by random walkers, 
which after being released at the periphery of the sys- 
tem, diffuse towards growing cluster and freeze on it. To 
simulate DLA, several numerical techniques have been 
developed, of which most powerful are the off-lattice al- 
gorithms 

Applications to other Laplacian problems have been 
proposed based on the random walks idea. In partic- 
ular, handling problems such as viscous fingering within 
this framework requires reducing noise of individual walk- 
ers as well as modeling surface tension. Reduction of 
noise was achieved by the method of multiple hits [^]|]] 
in which particles freeze on a particular site adjacent to 
the already grown cluster only after this site has been vis- 
ited more than n m i n times, where n m i n is an acceptance 
threshold. The method of noise reduction MM has been 
introduced in the context of the on-lattice DLA mod- 
els. More recently, this method was combined with the 
off-lattice technique and studied theoretically fl4] , |l^ , |l8f . 
In an advanced version p4[ of the multiple hits method 
the random walkers move off-lattice and sticking rules 
are defined by using a finite number m of antennas at- 
tached to each particle, where, for instance, m = 4 for the 
square lattice DLA. Each of m antennas has a counter 
which scores the number of times n random walkers ar- 
rive on it and is then used to set a threshold n m j n for 
freezing. Having a finite number of antennas used per 



each particle makes the sticking rules in the mutiple 
hits models anisotropic. This anisotropy is essential be- 
cause it is significantly amplified by the growth dynam- 
ics in the low noise regime of n m i n ^> 1. This built-in 
anisotropy of growth rules has been used to test univer- 
sality of DLA and also to simulate dendritic crystal 
growth @Q. 

It has been proposed that surface tension can be mod- 
eled by introducing a probability t < 1 of freezing upon 
each encounter with the cluster |lOj or by making freezing 
dependent on the local neighbors configuration jll]||[12| . 
In the model Jl(j] with t -C 1 each randomly walking 
particle freezes only after encountering the cluster about 
t^ 1 3> 1 times. As a result, the freezing point is displaced 
from the point of first encounter by the distance d ~ i" 1 
in the units of particle size. Effectively, in this model 
a finite length scale d is introduced over which the har- 
monic (Laplacian) measure describing the probability of 
the first encounter is being probabilistically averaged. It 
has been conjectured in the works |0,^,E2| (and partially 
confirmed by various features observed in the growth pat- 
terns) that the length scale d simulates capillary radius 
in the Laplacian problem with surface tension. 

The field of Laplacian growth, despite being well devel- 
oped by now, contains several long-standing unresolved 
problems. Firstly, the lattice simulation, although ex- 
tremely efficient algorithmically, does not seem to be a 
natural starting point for analytic understanding of large 
scale phenomena, such as fingering, fractalization and 
scaling. Secondly, the methods of simulating Laplacian 
growth that have been used so far are not entirely free of 
problems, the most notable being an intrinsic anisotropy 
of growth rules. Already the original DLA rules |2| use 
square lattice and thus are anisotropic. This anisotropy 
is weak and reveals itself only in very large DLA clus- 
ters P,0. However, when the noise level is reduced using 
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multiple hits, the underlying lattice anisotropy is ampli- 
fied @,||. The noise-reduced growth remains vulnera- 
ble to anisotropy even when off-lattice random walks are 
used, due to the anisotropy of freezing rules mentioned 
above. 

It was demonstrated recently that Laplacian growth 
can be studied using an entirely different approach based 
on iterated conformal maps |[9| . The model uses an- 
alytic functions chosen in such a way that upon acting on 
a unit circle they produce bumps of prescribed size. Iter- 
ated n times with the parameter defining the bump size 
chosen according to a certain rule, these maps produce a 
cluster of n particles of nearly equal size. The conformal 
model of growth has become recently a subject of active 
work ||||. 

The goal of this article is to extend the conformal map- 
ping methodology to the problems with reduced noise 
and growth anisotropy. Here one clear advantage is that 
the growth rules using conformal mapping |19(| are in- 
trinsically isotropic. Because of that one can easily avoid 
problems pertinent to other models, in which growth 
anisotropy and reduced noise are intertwinned. The 
main idea behind the noise reduction method proposed 
in this work is to average the Laplacian measure over fi- 
nite length which is larger than the particle size in the 
original model [flQfl . For that we alter particle shape and 
use "flat" particles extended along the cluster bound- 
ary and thin in the growth direction. To compare our 
method of reducing noise to other techniques, we note 
that the positions of flat particles are chosen strictly ac- 
cording to Laplacian measure, like in the multiple hits 
method H|| . The control over noise is achieved by sup- 
pressing noise at the length scales shorter than the parti- 
cle larger dimension. This is in contrast with the multiple 
hits method, where noise is suppressed due to statistical 
averaging over many particle growth attempts uniformly 
over all length scales down to particle size. Because of 
the appearance of a new length scale our method some- 
what resembles the surface tension models used in the 
DLA lattice growth 0Jll|,|,|l|. 

One notable difference from previous models is in the 
dependence of the computation time on the achieved level 
of noise reduction. Reducing statistical fluctuations in 
the multiple hits models required increasing the number 
of random walkers used to grow the cluster inversely with 
the noise reduction parameter. In our method one can 
reduce noise arbitrarily without increasing computation 
length, simply by varying particle aspect ratio with par- 
ticle areas kept fixed. We also demonstrate that growth 
anisotropy can be naturally incorporated in the confor- 
mal mapping method without affecting noise reduction. 

Our plan in this article is as follows. We start with 
revisiting conformal mapping model. We discuss some 
issues ignored before and propose a generalization to the 
problems with reduced noise and anisotropy. In Sec. [l| we 
review the model, keeping focus on aspects that will be 
important in the rest of the article. In Sec. |ni] we study 
the distribution of particle areas produced by growth 



rules employed in Ref. p9[ . We observe that these rules 
lead to occasional appearance of very large particles. 
To fix this problem, we evaluate particle areas at each 
growth step and apply an acceptance criterion for newly 
grown particles according to their area. In Sec. IV we 
describe a model with reduced noise. To suppress noise 
we use particles which are thin in the growth direction 
and smooth at corners. In Sec. we show how these 
growth rules can be generalized for anisotropic growth. 
In Sec. VI we study scaling properties of all introduced 
varieties of the model and compare them with each other. 
We find that the fractal dimension estimated from cluster 
radius scaling is less sensitive to noise than to anisotropy. 
For isotropic growth, both with and without noise reduc- 
tion, the dimension is very close to 1.7. For anisotropic 
growth, reducing noise to the level at which anisotropy 
reveals itself strongly shifts fractal dimension to some- 
what lower values. In this regime, the fractal dimension 
depends on symmetry, and is found to be 1.62 for the 
four- fold symmetry and 1.5 for the three- fold symmetry. 
Finally, in two appendices we discuss in detail the parti- 
cles area distribution and possible improvements of our 
numerical procedure. 



II. CONFORMAL MAPPING MODEL 

We describe growing cluster by a sequence of domains 
T>q C T>i C T>2 C ... corresponding to subsequent growth 
steps in time. In the canonical formulation growth 
occurs due to particles diffusing from infinity one by one 
and freezing as soon as they reach the cluster boundary. 
The particles concentration u(r) obeys diffusion equa- 
tion, which in the quasi-stationary approximation of slow 
growth is written as 



V 2 u(r) = with u(r) 



0, 



Irl — > oo . 



(1) 



Zero boundary condition on the cluster D„_i describes 
freezing of the n-th particle upon arrival with proba- 
bility one. The points of the cluster boundary <92?„_i 
where subsequent additions are made are selected ran- 
domly with the probability given by the so-called har- 



monic measure 



dP=\Vu\dl, (2/ C <9£> n -i : 



(2) 



where dl is boundary element of the cluster T> n _\. As 
the domain changes, ... — > T> n — > T> n+ i — > the prob- 
lem ([!]) has to be solved again for every new domain to 
determine from (|^) the new particle position probability. 

A considerable computational simplification of the 
problem can be achieved []l9| by using a sequence of an- 
alytic functions F n (z), n = 0,1,2,..., to represent the 
domains T> n . The functions F n are chosen so that each 
of them defines a conformal one-to-one mapping of the 
unit disk | z\ < 1 on the domain T> n , including boundary. 
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Adding a new object to the cluster at the n-th growth 
step is described by changing the mapping F n as follows: 



Fn(*)=J r n-l(/A B ,fl B (*)), F {z) = Z. 



(3) 



Here the function fx n ,9„(z) maps the unit circle \z\ = 1 
onto a unit circle with a bump centered around the point 
z = e l6n of the circle. The bump size is determined by 
the parameter A„ as discussed below. The angle 9 n is 
chosen randomly at each growth step. 

The key simplification that arises in the conformal 
mapping representation (^) is due to the fact that the 
harmonic measure (|^) is translated into a uniform prob- 
ability distribution for 9 n , so that dP{9) = c19/2tt. Also, 
there is no statistical correlation between subsequent #'s. 

The form of the function fx,e(z) growing bumps can 
be chosen according to the computational needs |19| . In 
this article we use 



x.o 



(z) = e i6 g- 1 (fx(g(e- i6 z))) 



(4) 



where the function g{z) = [z — l)/(z + 1) maps the unit 
disk \z\ < 1 onto the left half- plane Rez < 0, and the 
function 



fx(z) = hx(z)/h x {l) , h x {z) = z + Vz^+X 2 (5) 

grows a semicircle of radius A, as shown in Fig. [I]. The 
function fx(z) is defined in (||) so that /a(1) = 1. This 
ensures that the mapping (|]) maps z — oo onto itself. 




FIG. 1. The sequence of mappings constituting 
f x ,e(z) defined by @, (§). 

Ideally, the values A„ defining particle size should be 
chosen so that all particle areas are equal. In the confor- 
mal mapping model ]l9| this is approximately realized 
via predicting the bump size to be obtained at the n- 
th step using the Jacobian of the already-grown cluster 
mapping F„_ x . 

The argument is as follows. The area of the semicircu- 
lar bump grown using fx n (z) is 7rA^(l + 0(A 2 ))/8. The 
area of the corresponding bump produced by F n (z), at 
small A„, is approximately | J n -i | 2 7rA 2 /2, where J n -i is 



the Jacobian of F„_i evaluated at the position of the 
n-th bump: 



Jn-l — F'n-\{ z 



(6) 



Hence, to compensate for stretching due to the Jacobian 
J„_i, one has to choose the values of A„ as follows: 



A n — | Jn-l | 1 Aq 



(7) 



where the parameter Ao defines particle size. For the 
growth involving particles of very small size the rule (0) 
would have been sufficient to ensure identical areas of 
all bumps. For our problem, in which the bump sizes 
are small but finite, the areas are only approximately 



equal. However, 



one can demonstrate that, 
the 



improvements discussed in Section III 



after certain 
rule (0) pro- 
duces bumps with sufficiently close areas. 

The form (Q) , (||) of the mapping / has several nice fea- 
tures. First, since fractional linear function g(z) maps a 
circle onto a circle, the mapping fx.e produces a crescent- 
shaped particle with circular boundary (see Fig. [I]) . Sec- 
ond, a simple calculation shows that the particle curva- 
ture radius equals A. The latter has the following conse- 
quence. Consider growth starting from a circular cluster 
of radius r, described by the mapping F(z) — rz. The 
mapping F(fx t g(z)) then produces a particle of curvature 
radius equal to Ar. After the value of A is chosen accord- 
ing to the rule (0), A = Xq/F' = Xa/r, the particle radius 
becomes equal to Ao, independently of the cluster radius 
r. The area of this particle is readily evaluated: 

a(X ,r) = a* + A r - (r 2 - A 2 ,) tan _1 (A /r) , 

a, = 7rAg/2. (8) 

The area a(Ao,r) varies between a* for Ao -C r and 2a* 
for Ao 3> r. For a generic non-circular cluster the particle 
area cannot be found analytically. Statistics of the areas 

We will see that typical 



will be discussed in Section III 



area of a particle is of the order of a* . 

The overall size of the cluster T> n grown according to 
(H); (0) i s wen characterized |Q by the mapping F n 
stretching factor at large scales: 



Rr, 



F' n {z oo) = JJ f'x k .e k ( z -> °°) 



(9) 



k=l 



The cluster radius R n can be conveniently evaluated us- 
ing (j^) together with the property 5(00) = 1 as follows: 



Rn 



n 



2 = 1 



2\l/2 



(10) 



k=l 



k=l 



The reason for R n to be an accurate measure of the clus- 
ter T> n radius lies in the properties of so-called univalent 
functions @|l|. 

At large n the cluster radius R n is expected to grow as 
n Q , where a is a numerical constant. This is consistent 
with (ffi| provided that A 2 n -> 2a at large n HEifJ]. The 
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growth problem (Q), (0) is believed to give rise to fractal 
objects with fractal dimension d < 2. There are several 
conventional definitions of fractal dimension of a growing 
cluster JjJ] . In this article we employ scaling of the cluster 
size with its area. Also, one can use box counting, or the 
relation between average growth velocity in a strip and 
the strip width. 

Taking R n defined in (|^) as a cluster radius provides 
a numerically efficient method for calculating fractal di- 
mension. For that one looks for a scaling relation of the 
form 



Rn « A)( d , 



(11) 



where A n is total area of the cluster T> n . The dimension d 
is related with the parameter a describing scaling of A„ as 
follows: d — a -1 , which is true provided A n oc n. In our 
simulation we make sure that individual particle areas 
have a sufficiently narrow distribution (see Section III), 
and thus total cluster area is indeed proportional to the 
particle number. 

Scaling properties of the growth problem described 
above have been explored by several groups ]l9|-p5[. It 
was concluded that the properties of the growth result- 
ing from the conformal mapping model match those of 
the lattice DLA models. Below we revisit the relation 
between the problem (fi]), (||) and the conformal map- 
ping model (||), (||), ([?]) and discuss several interesting 
extensions of this model. 



III. AREA DISTRIBUTION 



In the original work |19| it was assumed that the rule 
(0) is sufficient to produce particles with nearly equal 
areas. This assumption was apparently consistent with 
the cluster images in which each particle is represented 
by one or few points. To investigate this issue more 
closely, in this work we have chosen a different method 
of representing particles, in which the exact boundary 
of each particle is shown. An example of a cluster with 
the boundaries of all individual particles displayed (see 
Fig. ||), demonstrates that the areas of almost all of the 
particles are indeed quite close. However, there is also a 
number of exceptional particles of large areas. 

Large particles tend to appear within fjords and seal 
the space between well developed branches. Typically, 
this happens when particle growth is attempted on the 
periphery of an actively growing region. Insufficiency of 
the rule (Q) for keeping particle areas small is caused 
by fluctuations of the Jacobian F^(z) over the unit circle 
\z\ = 1. These fluctuations can be large in magnitude and 
also very abrupt, happening on the scale of the order of 
X n within the circle arc mapped onto the particle bound- 
ary. In the case when a newly grown particle overlaps 
with such a fluctuation, it can be "artificially stretched" 
under the mapping. 



The appearance of large particles has been reported 
in Ref. |2lj] and a method for eliminating them was pro- 
posed, based on choosing an optimal shape of particles 
produced by the mapping f\ t g(z). It was argued |^T| that 
the best value of the parameter < a < 1 in the map- 
ping defined in Ref. |lij] is given by a — 2/3. This value 
provides a compromise between abundance of large par- 
ticles at a — > and needle-like particles shape at a — > 1. 
Since in this article the particle shape will be used to 
tune noise, we employ a different method for eliminating 
large particles, as described below. 




FIG. 2. Cluster of N = 400 particles grown using the 
model ||, (@) with Ao = 0.2 and f x ,e of the form (§), 
(^). The boundary of each particle is displayed. Note 
large particles which appear rarely and seal fjords. 

To study the role of exceptional particles, we calculate 
the particle areas generated by the growth model (||) , (pi). 
The method employed to evaluate particle areas is the fol- 
lowing. For the particle grown at a step n, we start with 
few points on the unit circle \z\ = 1 which are mapped by 
F n on the particle boundary. Subsequently, we add new 
points on the circle \z\ = 1 in between the old points, and 
compute distances between images of neighboring points 
under all mappings 



fn 7 fn — 1 ° fn ) 



F, 



/i°/ 2 o...o/ n , (12) 



where is a shorthand notation for f\ ki e k and ° stands 
for mapping composition. We keep adding new points 
until the distances between the images of all neighbors 
will not exceed 7A0, where 7 <C 1 is a numerical fac- 
tor. We used the above procedure with 7 = 0.05, which 
produces about two to four hundred points per particle. 

This method enables one to have an accurate graph- 
ical representation of each particle, as demonstrated in 
Fig. ||, and also to evaluate particle areas with the accu- 
racy on the level of 0.1%. This is illustrated in Fig. g(a) 
showing how the cluster total area A n is changing dur- 
ing the growth of the cluster displayed in Fig. ||. The 
area A n grows as a function of n in small steps of or- 
der a* = 7tAq/2, alternated with occasional jumps of a 
much larger magnitude. The jumps correspond to the 
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appearance of large particles which seal inner cluster re- 
gions. The decomposition of the growth of A n into the 
smooth and singular parts is revealed more clearly in 
Fig. ||(b) showing the dependence of A n versus n for the 
same growth as in Fig. ^(a), with a linear function 2.1 a* n 
subtracted. 




(b) 




400 



200 



400 



total cluster area A n dy- 
namics for the growth of the cluster in Fig. |2| 
N = 400, Ao = 0.2; Dotted line: the area 
An versus n for the growth in which only parti- 
cles with the area inside the window [0, 3a,] are 
accepted, other growth parameters are the same; 

(b) The same area dynamics as in (a) with a linear 
function subtracted: A n — 2.1 a,n versus n. 

A histogram of the individual particles areas a n = 
A n — A n -i is plotted in Fig. ^(a). The area distribution 
was obtained by averaging over 10 realizations of the first 
1000 growth steps with the parameter Ao = 0.2 (same as 
in Figs. [| ||). The area distribution V{a n ) is sharply 
peaked about 2.1 a* and has a number of other interest- 
ing properties that will be discussed in Appendix |A|. 

The feature of the distribution V(a n ) displayed in 
Fig. [|(a) which corresponds to the exceptionally large 
particles is the tail stretching far to the right from the 
main peak. Note that only the beginning of the tail is 
displayed in Fig. [|(a) because the weight of the tail in the 
probability distribution is insignificant, and so the values 
of V(a n ) far in the tail are too small to be visible on the 
scale of the peak. 

To display the tail we replot the distribution V(a) on 
a log-log scale, as shown in Fig |. The right tail of V(a) 
is power-like, V(a) oc a~ M with fi w 2.5. Since /i > 2, 
the first moment (a) = JaV(a)da is finite (and thus 
(a) ~ a*). However, since [i < 3, the second moment 
of Via) is divergent. The existence of the mean particle 
area (a) means that (A n ) oc n. However, the absence of 
variance implies that the fluctuations of A n about the 
mean value are non-Gaussian and larger than required 
by the central limit theorem. Both features are clear in 
the sample A n dependence in Fig. ||. 

Let us remark that the tail in Fig. || in its far end is 
apparently somewhat steeper than a _M . We believe that 
this deviation from the a _M behavior is due to finite num- 



ber N — 2000 of time steps in the growth samples used 
to obtain V(a). In a finite cluster there is an upper cut- 
off on possible particle areas. This makes the far tail of 
V{a) non-stationary, shifting the cutoff to larger areas as 
N increases. 





a/a* 



a/a* 



FIG. 4. (a) Probability distribution function of par- 
ticle areas normalized by a, = 7rAo/2 ~ 0.063 (see 
(§))• Statistics was taken over 10 independent runs of 
the growth with N = 1000, A = 0.2. Note the peak 
at ~ 2.1a, and the tail corresponding to exceptionally 
large particles. The largest area observed was « 87.8 a* 
(see Fig. ^|). Small peak marked by arrow is due to 
the primary particles growing directly on the unit circle. 

(b) Same as in (a) for the growth in which only parti- 
cles with areas in the window [0, 3a,] are accepted. 
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In a/a* 

FIG. 5. Logarithm of the probability distribution 
function of particle area In V(a/a t ). Calculated from 
83 realizations with N = 2000, A = 0.8. 

Clearly, one would like to inhibit the appearance of 
large particles with areas in the tail of V{d). This is 
desirable because, even though the tail is quite thin and 
large particles are rare, occasionally appearing extremely 
large particles may affect macroscopic characteristics of 
the growth. In particular, the relatively slow power law 
decrease in the tail, V(a) oc a~ M , may affect scaling of 
the cluster size R n and/or the numerical accuracy of the 
scaling exponent. To eliminate the growth of large par- 
ticles, we choose an acceptance threshold a max = 3a* to 
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truncate the tail of the distribution in Fig. |^(a). Then, 
for each growth step, we calculate the new particle area 
a n . The particle is accepted only if a n < a max , otherwise 
the particle is discarded and a new attempt of particle 
growth is made. An example of the cluster grown accord- 
ing to these rules is displayed in Fig. ^]. 

One can notice immediately that the overall structure 
of the branches in Fig. ^| is much more regular than that 
in Fig. ^. The distribution of areas for such a growth 
is shown in Fig. |^(b). Within the acceptance window 
[0, flmax] the distribution V(a) repeats in all details the 
distribution shown in Fig. 0(a) for the growth with all 
particle areas accepted. 




FIG. 6. Cluster grown with the particle area ac- 
ceptance window [0,3a*] (see Fig. §(b)), N = 17545, 
A = 0.8. 



For the growth with large particles eliminated, the area 
A n increases as a linear function of the step number. Av- 
erage increment of A n is given by the mean value of a n 
taken from the distribution shown in Fig. ^(b) . To verify 
this, we plot A n versus n in Fig. for the same growth 
parameters as those used in Fig. g7 where only the areas 
up to 3a* are accepted. Note a small difference between 
the slope of the dependence at n < 50 and at larger n 
that appears because of relatively smaller size of the pri- 
mary particles growing directly on the unit circle. 

The growth model augmented with the area accep- 
tance criterion has a new parameter a max /a». In prin- 
ciple, choosing different values of a max gives rise to dif- 
ferent growth patterns. However, as long as the window 
[0, a m ax] contains much of the V(a) peak area, we do not 
observe any qualitative change in the growth. 

Scaling properties of the growth can be studied in 
several ways. Previous studies of scaling [l9|j2l| -p4]| are 
based on the relation R n cx n a , where R n was obtained 
for the growth with unrestricted areas. However, it would 
be more in the spirit of the notion of a fractal to use 
the relation (|ll]) between cluster size and its area. This 
clearly would not work well with large particles being 
present, because statistical fluctuations of the cluster 
area A n are quite large in this case. On the other hand, 
for the growth with restricted areas used in this work, 
the fluctuations of A n arc reduced to the level consistent 



with the central limit theorem, and thus one can employ 
the relation ( [Tl"|) to study scaling. 

Of course, it is not clear a priori whether the growth 
with area cutoff is equivalent macroscopically to the 
growth with unrestricted areas. On general grounds, one 
may expect the growth to be significantly altered by elim- 
inating large particles. Whether this is true can be in- 
directly tested by comparing the R„ vs. n dependence 
for the growth with unrestricted areas with the fractal 
dimension obtained from the relation ([ll]) for the growth 
with area cutoff. We note in that regard that the scaling 
exponent d = 1.7 found below (see Section fVl| ) matches 
exactly the value found in Ref. [^9[ . However, although 
the presence or absence of large particles seems to be 
irrelevant for the cluster size scaling, other growth char- 
acteristics, such as the structure of branches and fjords, 
are likely to be more sensitive to the method of treating 
large particles. 

We postpone the discussion of various details and fea- 
tures of the area distribution "P(a) to Appendix [A| In the 
remaining part of the article we use the conformal map- 
ping model augmented with the area acceptance criterion 
to study several interesting Laplacian growth problems. 



IV. NOISE-REDUCED LAPLACIAN GROWTH 

Roughness of the growing cluster is mainly due to two 
factors: shot noise and the Mullins-Sekerka instability 
p7| . The shot noise results from the randomness of the 
aggregating particles positions, and so it contributes to 
the fluctuations equally on all spatial scales down to the 
particle size. The Mullins-Sekerka instability is due to 
aggregation rate enhancement near the tips, which leads 
to incremental growth of perturbations of a smooth front. 
The wavenumber dependence of the growth rate for a har- 
monic modulation of an interface moving with average 
velocity v is given by 7^ = v\k\. The linear fc-dependence 
of 7fc implies that the instability develops first on the 
smallest scale, in our problem given by the particle size. 

To study the ultraviolet cutoff role, i.e., the effect of 
short distances on the noise and the instability, it is of 
interest to introduce a parameter in the problem which 
allows to shift the value of the cutoff scale to values larger 
than the particle size. One expects that upon doing so 
both the noise and the instability growth rate will be 
reduced. 

In the mapping model, the noise level can be controlled 
by altering the shape of aggregating particles. Below we 
show how by changing the function f\(z) 7 defined by (|^), 
one can create "flat" particles which are wide along the 
interface and thin in the growth direction. The reason 
that noise is suppressed due to using flat particles is the 
following. In this growth, a particular displacement of 
the growing cluster boundary amounts to a larger num- 
ber of layers than in the case of rounded particles used in 
Ref. p9[| . Then, due to statistical averaging over many 



G 



particle layers the boundary displacement becomes less 
erratic, and so the noise is reduced. Quantitatively, the 
noise suppression factor can be estimated as a square root 
of the particle aspect ratio. 

Flat particles can be produced by modifying f\(z) as 
follows: 

h. p (z) = Whil OMMz))) , \ = z~r~nz > (13) 



W = h 



p > 1. 



(14) 



The function h\(z) is defined in (|^), and its inverse has 
the form h' x l {z) = \{z — X 2 /z). The factor W is in- 
troduced in order to have /a, p (1) = 1, like for the func- 
tion f\ (z) defined by (||) above. The resulting function 
(^) satisfies f\^e(oo) — oo, which ensures the property 
F n (oo) = oo for all n. 

The mapping produced by the function (O) is illus- 
trated in Fig. 0. Note that because of h\ p (i\)/p = i\ p , 

the square root singularities in f\, p at z — ±i\ are absent 
for allp > 1. Instead, the mapping composition ( |l3| ) pro- 
duces weaker singularities of the form (z ± iX) 3 ^ This 
smoothens the corners of the particles, as shown in Fig. |t]. 




^ FIG. 7. The sequence of mappings constituting 
f\, p (z) as defined by @, p = 1.5. 

Qualitatively, under variation of p the particle shape 
evolves as follows. At p = 1 the mapping f\ p form Jl3| ) 
coincides with (|^). Increasing p produces particles with 
growing aspect ratio, as can be seen from comparing the 
zoom parts of Figs. |To[ [ll] and ||. 

To illustrate the effect of p on the particle shape, con- 
sider the mapping function (|13| ) in the limit p> 1. First, 
one can rewrite (fl3) as 



A,p(«) = W h x (z) 



A> 2 - 1) 



(15) 



where W' is a prefactor chosen so that /a,p(1) = 1- Ex- 
panding ( IHf ) to lowest order in 1/p, one obtains 



f\ P (z) 



z H r 



A 2 / z + 2Vz 2 + A 2 



(z + \fz r 



A 2 ) 2 



Cz 



(16) 



where C* = (2h A (l) - l)//i 2 (l) and /i A is defined in (|). 
The boundary of the particle produced by f\ jP (z) of the 



form ([T^), to lowest order in 1/p, is 



A 2 p" 



(A 2 



y ) 



(17) 



where x + iy = z. The area of this particle is 37rA 2 /4p 2 . 
Mapped by g^ 1 , according to the area is multiplied 
by a factor equal to 4 at A -C 1. 

One can use the growth mapping model (|§|), (Q) with 
the new function ( |13| ) to grow clusters in a pretty much 
the same way as it was done for the model with p = 1 
in Section III. The first step is to study the particle ar- 



eas distribution for the growth with unrestricted areas. 
The distribution looks similar to that in Fig. |], contain- 
ing central peak and the tails corresponding to very large 
and very small particles. In this case the peak is some- 
what wider than for the p = 1 case. However, much of its 
weight in the distribution V(a) is still contained in the 
window [0,3a*]. Here the "standard area" a* is defined, 
by analogy with (||) , as the area of a particle grown over 
a perfectly flat interface. (For p ^ 1 there is no closed 
form expression for the particle area, like (||), and so one 
has to calculate a* numerically.) 

As before, at each growth step we choose 6 n randomly, 
< n < 27r, and calculate the parameter A„ using (|7|), 
i.e., based on the particle area predicted from the Jaco- 
bian J n -i- Then we evaluate the actual of the 

particle. To inhibit the appearance of large particles, we 
use the acceptance window [0, 3a*]. If a n > 3a*, the par- 
ticle is not accepted and a new growth attempt is made. 




FIG. 8. Cluster grown with p — 
number of growth steps N = 15043. 



3, Ao 



0.2. The 



An example of growth with p = 3 and Ao = 0.2 is 
displayed in Fig. ||. In the inset we zoom on the details 
of one finger. Note that individual particles are indeed 
quite flat and are evenly spread over the cluster bound- 
ary, indicating reduced noise. The growing interface is 
overall very smooth, without sharp tips or corners. Also, 
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the fingers are much thicker than for the p — 1 growth 
(see. Fig. ||). 

The cluster size R n is defined by (51) . As in Section II] , 
the terms in the product (|J) can be evaluated using the 
relation f' Xp ^(oo) = l/j' Xp (l), where 

? h x (i)/VT+Vhl p (h x (i))+ p 2\l 



^/j2(l) + A2 ^a p (^a(1))-P 



2 A 2 



In the following Section VI we use @ along with (|) to 
evaluate the cluster radius R n and study its scaling. 

The appearance of the cluster in Fig. || shows that us- 
ing flat particles indeed helps to reduce statistical fluctua- 
tions. In this model, effective averaging of harmonic mea- 
sure is due to the presence of a tangential-to-boundary 
length scale set by particle larger dimension. This length 
scale is controlled by the parameter p and becomes large 
at p > 1, if measured in the units of particle size -yjal. 
Noise reduction takes place due to the absence of fluctu- 
ations with wavelength smaller than the particle larger 
dimension, resulting in the shift of the shot noise spec- 
trum cutoff wavenumber from 2ir/^/a^ to lower values as 
the parameter p is increased. Because of reduced noise, 
as compared to the p = 1 case, more agregation events 
of flat particles are needed to reach a given radius of the 
cluster. 

Averaging over a tangential length scale is somewhat 
similar to that used in the on-lattice DLA models to sim- 
ulate surface tension Jl(],|ll],[3 12 . In these works freezing 
of random walkers upon each ecounter with the cluster 
was described by a finite probability t < 1 which could 
be a function of occupancy of the sites around freezing 
point. Since freezing of each particle typically takes place 
after about t~ x encounters with the cluster, at t -C 1 
these models are characterized by a large length scale 
over which Laplacian measure is probabilistically aver- 
aged. Similarly, the flat particles used in our model can 
be thought of as a result of averaging over possible par- 
ticle positions within a finite length scale taken over har- 
monic measure. Moreover, there is a slight dependence 
of particle size on their growth position: the particles 
appearing near the tips are somewhat smaller than those 
appearing in the concave regions (see Fig. |^). This cor- 
relation is consistent with the surface tension interpreta- 
tion. 

The crucial difference, however, is that particle posi- 
tions in our model are chosen according to the unaltered 
harmonic measure, whereas in the surface tension mod- 
els particle freezing depends on local boundary geometry. 
From that point of view our model is more similar to the 
multiple hits models [^],[)| in which statistical averaging 
of the harmonic measure over particle growth attempts 
is used to control noise. In these models noise reduction 
is achieved by averaging over independent random walk- 
ers with a threshold on the minimal number of visits of 
each site required before freezing at this site. Since inde- 
pendent walkers arrive at very distant points of cluster 



boundary, this averaging is not characterized by an addi- 
tional large length scale and thus bears no resemblance 
to surface tension. 

The models using finite freezing probability t < 1 have 
been shown to give rise to clusters with thick branches. 
The Laplacian character of the dynamics and the analogy 
of the averaging length scale with capillary radius was 
pointed out p0[ and a relation with the Saffman- Ta ylo r 
problem with surface tension has been conjectured p2| . 
Because of the large length scale appearing in our averag- 
ing scheme, here a similar relation to the problems with 
surface tension can be conjectured. Indeed, the growth 
displayed in Fig. || looks like a typical fingering pattern 
observed in the Saffman- Taylor problem with surface ten- 
sion. As a word of caution, one should realize that all 
available evidence for the equivalence between the prob- 
lem with surface tension and our large p growth, however 
similar they appear to be, is rather indirect. The issue 
of whether or not this growth model is indeed character- 
ized by an effective surface tension and how the latter 
compares to the noise will be discussed elsewhere. 



V. ANISOTROPIC GROWTH MODEL 

The iterated mapping model (||) , ([7j) can be generalized 
to describe spatial anisotropy of local growth rate. Such 
an anisotropy is characteristic of crystal growth, in which 
all particles arriving at the crystal-liquid interface have 
to accomodate to the anisotropic crystal structure Q. 

Anisotropic growth oftenly gives rise to anisotropic ir- 
regular fingering patterns called dendrites 11 ll| |l7| , |i~| . 
The dynamics of dendrites growth obeys scaling laws sim- 
ilar to that of the isotropic Laplacian growth 0,^8| . One 
of the outstanding theoretical questions is how the scal- 
ing exponents depend on the anisotropy. 

In this problem, the cluster grows due to spatially 
isotropic diffusion and aggregation of particles. Thus 
the quasistationary probability distribution still obeys 
Eq. (|lj). The new element, compared to the isotropic 
model, is that the cluster change due to particle freezing 
at the boundary depends on the local growth direction 
v, |v| — 1. (The unit vector v is normal to the bound- 
ary.) Accordingly, the probability of joining the cluster 
becomes a function of v, and Eq. (0) is replaced by 



dP = fi(v) \Vu\dl, dlcV„ 



(19) 



where the function il(v) describes anisotropy. 

In order to include anisotropy in the mapping model 
(H) i (@) , at the n-th growth step one has to be able to pre- 
dict local growth direction v„ from particle positions de- 
scribed by randomly chosen angles 9k, k = 1,2, ,..,n— 1. 
This is possible because complex-valued Jacobian of a 
conformal mapping keeps track of the angle change under 
the mapping. Specifically, consider 9„ = 9 n + arg J n _i, 
where J„_i is given by (^). Then 0„ defines a normal 
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to the cluster boundary, v n = cos Q n x + sin <d n y, at the 
growth point F„_i(e "). 

Now, there are several possible ways to account for 
the growth anisotropy. For instance, one can introduce 
the anisotropy by making A n a function of 0„, e.g., 
A„ oc ri 1 / 2 (v„). Another way is to introduce acceptance 
probability for the particles which depends on O ra in some 
way. In the simulations reported below we use an accep- 
tance window for MQ n with M = 3,4, corresponding 
to the growth with an M-fold symmetry Only particles 
with ©„ such that 



< M6„ < 



(20) 



are accepted. Here 9 max is a parameter in the interval 
[0, 7r] controlling the degree of anisotropy. Small val- 
ues # m ax <C 7r correspond to highly anisotropic growth, 
whereas fully isotropic growth is recovered in the limit 

Other aspects of the simulation are the same as in Sec- 
tion IV. Wc employed the elementary mapping f\ tP (z) 
of the form ( p"3| ) with the noise level controlled by the 
parameter p > 1. Particles with large areas were elim- 
inated using the acceptance window [0, 3a*] defined in 
Section IV. An example of growth with the three- fold 
symmetry (M = 3) is shown in Fig. ^. In this case, we 
used A = 0.8, p = 1.5, and max = cos- x (0.9) « 0.451. 
The cluster is characterized by overall symmetric main 
branches covered with numerous side branches. 

In our model one has a separate control over the degree 
of anisotropy and over noise, via the parameters 9 max and 
p. This is convenient for studying the effects of noise on 
the ordering of branches in dendrites. To illustrate that, 
we compare t wo g rowths with four-fold symmetry, dis- 
played in Figs. [l(J |ll]. The cluster in Fig. [l(] is obtained 
using the p — 1 model without noise reduction, as de- 
scribed in Sections EL H. The cluster in Fig. ftljis grown 
using the noise reduced model of Section IV with the pa- 



rameter p = 2. In both cases, we use the same anisotropy 
parameter: 6 max = cos -1 (0.95) ~ 0.318. One notes high 
anisotropy of the growth present at small scales in both 
cases, which is significantly suppressed at larger scales 
for the noisy growth with p = 1 (see Fig. [To| ). However, 
the p = 2 growth with low noise remains very anisotropic 
at all scales (see Fig. |TT| ) . 

It is known from studies of on-lattice DLA models that 
noise, no matter how strong, gives way to anisotropy at 
sufficiently large scales [p|p| JT7Pl5| ] . We thus expect that 
a similar effect may take place in the noisy growth with 
p = 1, making the growth shown in Fig. |l0| at larger 
scales looking like that in Fig. [n} 

In agreement with the studies of the off-lattice DLA 
models yl|, we observed that the dendritic growth with 
the symmetry of order M = 3, 4 is much more stable 
with respect to noise than that with M = 5, 6 or higher. 
Scaling properties of the anisotropic growth will be stud- 
ied in Section VI. 





FIG. 9. Anisotropic growth with M = 3 obtained 
using the window ( p(i| ) for the growth direction with 
cos # m ax = 0.9. Other parameters used: iV = 10146, 
A = 0.8, p = 1.5. 




FIG. 10. Anisotropic growth with M = 4, 
cos 6 max = 0.95, and high noise: p = 1, N = 7635, 
A = 0.3. 




FIG. 11. Anisotropic growth with M = 4, 
cos # max = 0.95, and low noise: p = 2, N — 6782, 
A = 0.8. 
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VI. SCALING PROPERTIES 



Scaling of R n for all growth models introduced above 
is studied here using the following procedure. The clus- 
ter radius i?„ obtained from (^), (]T(i|), (18), is plotted 
against the cluster area A n , evaluated as a sum of in- 
dividual particle areas a n . Asymptotically, at large n, 
one has Rn oc A\[ d . To determine d more accurately 
we optimize initial conditions of the growth, represented 
in our model by nondimensionalized particle size Ao, as 
described below. 

In the log-log plot of R n versus A n one can clearly 
distinguish two regimes, initial growth and developed, 
or regular growth, characterized by somewhat different 
slopes of the corresponding parts of the lni? vs. \nA 
curves. Geometrical meaning of these regimes is as fol- 
lows. For isotropic growth, with or without noise sup- 
pression, the cluster initially consists of branches growing 
practically independently. Later, at the regular growth 
stage, the number of main branches is reduced to four 
or five, all interacting and competing with each other. 
For anisotropic growth with the M-fold symmetry the 
number of main branches is M at all stages of growth. 
Regular growth in this case is distinguished by many fin- 
gers appearing on the sides of M main branches. 

The initial stage is more pronounced when the particle 
size, determined by the value of Ao , is much smaller than 
the unit circle from which the growth starts. Since we 
are interested in the regular growth scaling, in each case 
studied we tried to optimize the value of Ao to shorten the 
initial growth stage, carefully checking that the variation 
of Ao has no detectable effect on the asymptotic slope 
of the lni? vs. In A curve. The benefit of shortening the 
initial growth stage is that, at constant number of parti- 
cles, it leads to longer regular growth and thus allows to 
extract the scaling exponent with higher precision. The 
resulting curves are presented in Fig. [l2[ as described in 
the figure caption and below. 

The optimal value of Ao determined for the isotropic 
growth with p = 1 is close to Ao = 0.8. For the scal- 
ing analysis we used the growth displayed in Fig. [| in 
which N = 17545, A = 0.8, p = 1. In Fig. ||, it cor- 
responds to the lowest of the curves marked by a,±. To 
eliminate the effect of fluctuations at the initial stage of 
the growth, we also generated the curves a2, a3, and &4, 
by averaging lni? over 5, 10, and 50 growth realizations 
with N — 5000, 1000, and 200 time steps, respectively. 

For isotropic growth with reduced noise we analyzed 
two growths with p = 3: 

curve bi with N = 11611, A = 0.1; 
curve b 2 with N = 15043, A = 0.2. 
The curve b 2 corresponds to the growth displayed in 
Fig. H At low noise, the fluctuations of R n are quite 
small, which makes additional averaging over realizations 
unnecessary. 

Note that for the isotropic growth models the strategy 
of optimizing Aq works quite well, allowing to almost en- 



tirely eliminate the initial growth region. The scaling di- 
mension found from the slope of the best straight line fits 
is close to 1.7. To study the deviation from 1.7, we sub- 
tract from all curves the linear function lni? = In A/1.1 
and plot the result in the lower part of Fig. [l^. Note that 
upon this subtraction the curves for isotropic growth, 
with or without noise suppression, become nearly per- 
fectly horizontal. Estimate of the deviation from the best 
horizontal line fit shows that the value 1.7 is accurate 
within 1%. 
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FIG. 12. Log-log plots of R p /(a) (a) and R/A"> (b) 
vs. normalized area A /{a) for several clusters described 
in the text. Here f3 = 1.7, 7 = 1/1.7 = 1/(3, and (a) 
is the particle area averaged over the cluster. The plots 
(a) and (b) are connected by an affine transformation. 

For anisotropic models, we consider three different 
growths: 

curve ci with N = 10146, A = 0.8, p = 1.5, M = 3; 
curve c 2 with N = 7635, A = 0.3, p = 1, M = 4; 
curve c 3 with N = 6782, A = 0.8, p = 2, M = 4. 

These curves correspond to the growths displayed in 
Figs. [| [l(], and [ll], respectively. As above, we subtract 
the linear function lni? = lnvl/1.7. However, after this 
subtraction, the curves Ci and C3 retain some residual 
slope. Estimating it, we conclude that the best value for 
the fractal dimension is d w 1.5 for ci, and d w 1.62 
for C3. The latter value agrees with the values d » 1.58 
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and d w 1.63 of the growth with M — 4 reported in 
Refs. ||§. 

For the curve C2 corresponding to anisotropic growth 
with noise, after subtracting lni? — In A/1.7 we do not 
find any significant residual slope. It is possible, how- 
ever, that the dimension 1.7 corresponds to the crossover 
regime and changes to a lower value at larger TV. Simi- 
lar behavior is known to take place in the on-lattice DLA 
growth [Q, where the dimension 1.7 observed for not very 
large clusters crosses over to 1.63 at N ~ 4 • 10 6 . 

To understand possible sources of errors in determin- 
ing the fractal dimension from lni? vs. hi A curves, here 
we consider how R n and A n fluctuate with n. The fluc- 
tuations of \n.R n gradually decrease with increa sing n, 
as can be clearly seen in the lower panel of Fig. [lj. A 
convenient way to analyze fluctuations is to plot pairs 
(In Rn, In A n ) for particular N, repeating growth many 
times. In Fig. 13 we present results for 10 3 growth sam- 
ples and several values of N. The resulting clouds be- 
come more compact as N increases, indicating that the 
fluctuations of In R n and InAjv are decreasing. 



InR 



In Rn = ^ m /fe( z ~* 00 ) 



(21) 



fe=i 




In A 



FIG. 13. Clouds of points (InR, In A) corresponding 
to 1000 realizations for N = 25, 50, 100, and 200, 
A =0.8, p=l. 

Let us first discuss fluctuations of lnA n . The to- 
tal area A n is the sum of individual particles areas afc, 
k = 1,2, ...,n. Assuming that the areas a,k are indepen- 
dent or, more precisely, have only short correlations, one 
obtains a Gaussian distribution with the variance cx n. 
(As we argue below, there exist long negative correlation 
of particle areas, which may further reduce fluctuations 
of A n .) The fluctuations of In A n are simply given by rel- 
ative fluctuations 5A n /A n , which means that for large n 
the distribution of lnAn is also Gaussian, with the vari- 
ance oc n . 

On the other hand, the radius R n is a product ^) of 
stretching factors J^ 00 ^ = f' k [z — > 00). Since 
for all fc, the quantity R n grows monotonically, so that 
R n oc Al/ d at large n. Thus the noise due to fluctuations 
of jj, is °f a multiplicative nature. One can write 



which suggests that the distribution of lni?„ is Gaus- 
sian, i.e., the distribution of R n is log- normal. Indeed, 
the log-normal fit perfectly describes the statistics of 
R n , as demonstrated in Fig. |l4|. However, attempting 
a Gaussian fit produces an asymmetric distribution de- 
viating from the observed distribution of lni?„. Thus, 
even though the relative fluctuations of R n are small, 
the statistics is best described as log-normal. 




In J? 



FIG. 14. Probability distribution of In R for N = 100 
and N = 200 calculated from 10368 and 2862 growth 
realizations respectively, Ao = 0.8, p = 1. 

Naively, Eq. ( |2l| ) implies growth of the variance of 
lni? n with n. However, Figs. [II], [li] demonstrate that, on 
the contrary, the width of In R n distribution is decreas- 
ing with n. This nontrivial behavior was first mentioned, 
without explanation, in Ref. ||21fl - 

To rationalize the observed sharpening of the distribu- 
tion of R n , one can argue as follows. We note that the 
dynamics of R n is characterized by a negative feedback. 
Consider growth of a cluster which at the n-th step has 
a radius smaller than average. Then the Jacobian of F n 
is typically smaller than its mean value at this number 
of particles. In this case, according to (0), subsequently 
growing particles will have larger Afe's, and thus larger 
areas, until the cluster radius will approach the average 
value. The evolution of a cluster which at a certain step 
has a radius larger than average can be considered in a 
similar way. This long-time anticorrelation of Afe's sup- 
presses the fluctuations of R n . Also, it produces long 
negative correlation of particle areas. 



VII. SUMMARY 

To conclude, growth models using conformal mappings 
have large flexibility allowing for independent control 
over noise and growth anisotropy. We generalized the 
model (l^] by using flat particles to suppress noise. It 
is essential that these models lead to an intrinsically 
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isotropic growth with reduced noise, in contrast with 
other previously studied models. Also, we demostrated 
that favoring growth in certain directions can be used to 
simulate anisotropy of the growth rate. 

Having separate control on the noise and anisotropy, 
we have been able to analyze their effect on scaling prop- 
erties. We found that the fractal dimension d — 1.7, uni- 
versally for any isotropic growth, regardless of the noise 
level. However, the fractal dimension is somewhat re- 
duced in the presence of anisotropy. 

It was assumed [[l9| that particle size fluctuations, 
present in the conformal mapping model, are insignifi- 
cant. We observed that the growth rules used in Ref. [|Ts|] 
lead to occasional appearance of exceptionally large par- 
ticles. We have shown that by augmenting the model 
with an area acceptance criterion this problem is fixed. 

Clearly, more work has to be done to establish a re- 
lation of the introduced models with real physical pro- 
cesses, like viscous fingering or dendritic crystal growth. 
Another interesting open question is how to introduce an 
effective surface tension. 



ACKNOWLEDGMENTS 

We thank G.E. Falkovich, M.B. Hastings, V.A. Kaza- 
kov, V.V. Lebedev, B.Z. Spivak, and P.B. Wiegmann 
for useful discussions. Hospitality of Weizmann Institute 
made our collaboration possible. This work was partially 
supported by the Minerva Foundation (M.S.) and by the 
NSF Award 67436000IRG (L.L.). 



APPENDIX A: INSIDES OF THE PARTICLE 
AREA DISTRIBUTION 

Here we discuss in more detail the distribution of par- 
ticle areas. The main feature manifest in the area his- 
togram plotted in Fig. U is a sharp asymmetric peak at 
« 2.1a*. This peak has its origin in the dependence of 
particle size on the growth point. 

The argument is as follows. First we note that the 
growth is taking place predominantly at the tips of the 
branches. Because of that, for several particles growing 
on each other, there is a tendency to preserve growth di- 
rection. This leads to formation of relatively long chains 
of particles growing in a particular direction, clearly seen 
in the inset of Fig. ^. The chains are mostly formed at 
the tips of outer branches. 

Now, consider a particle growing near one of the tips. 
The area of this particle has some dependence on the po- 
sition of the growth point relative to the tip. The peak in 
the histogram in Fig. ^ is explained if one assumes that 
the particle area has a local minimum in the forward 
growth direction. The minimum in the area leads to a 
caustic in the histogram. Ideally, this would have pro- 
duced an asymmetric square root singularity with prob- 



ability equal zero on the left side. Because of particle size 
variation among branches, the singularity is smeared into 
a peak. 
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FIG. 15. Area of a particle as a function of its growth 
point, characterized by 6. 

To verify the above assumption, we consider areas for 
the first few particles grown on the \z\ = 1 circle with the 
parameter Ao = 0.2. The area of the very first particle is 
close to 1.2 a* and, according to (ph, is independent of its 
position. The area of the second particle a2 depends on 
its position 9 relative to the first particle, as shown by a 
solid line in Fig. |l5](a) . Note that the area is the same as 
that of the first particle when the particles are far apart, 
9 ^> A , and is overall substantially larger when the par- 
ticle overlap, 9 ~ Ao. Partially, this is explained by the 
dependence (||) of particle area on the circle radius. (As- 
suming that a 2 (9 ~ Ao) can be crudely estimated by (||) 
with r = Ao.) Another effect that contributes to the area 
ai increase for 9 ~ Ao is the variation of the Jacobian 
as a function of 9, leading to additional stretching of the 
second particle. 

The feature in Fig. |l5|(a) which is of interest in con- 
nection to the peak in the area distribution Via) is the 
minimum of a 2 (9) at 9 = 0. Translated to the histogram 
of areas, it leads to a caustic described by a square root 
singularity. However, as a possible explanation of the 
peak in Fig. ^ this is only partially satisfying, since one 
has to understand why similar caustics due to the two 
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maxima of a-i (9) are not observed in Fig. |j. 

The reason for the difference between the effects of 
maxima and minima can be seen from a comparison with 
the case of three and four particles. Consider the situa- 
tion when the second particle is centered exactly on the 
first particle, and the third particle is grown at an an- 
gular position 9 relative to the first two particles. The 
area of the third particle 0,3(0) is plotted in Fig. |l5|(a) in 
a dashed line. Note that, since the curvature at the min- 
imum of 03(0) at 9 = is much smaller than for a 2 (0), 
the corresponding caustic in V(a) will be much stronger. 
On the other hand, the curvature at the maxima of 03 (9) 
is about the same as that for 0,2(9). Both observations 
remain correct for any number of particles. To illustrate 
this we plot the area 04(8) of the fourth particle in the 
presence of three particles grown exactly on the top of 
each other — see dotted line in Fig. |l5|(a). 

Another notable feature in the plots of 02,3,4(0) is that 
the area becomes much smaller than a„ , approaching zero 
near certain values of 9. This behavior is related to the 
growth near particle corners, which are the points of di- 
vergence of the Jacobian. According to (0), larger Jaco- 
bian translates into smaller particle area. The particles 
growing near corners form the tail of the area distribu- 
tion V(a) at small areas a <C a*. The behavior of V(a) 
in this tail, V(a) cx a 1 / 2 , follows from the square root 
divergence of the Jacobian at particle corners. The slope 
1/2 is clearly seen in the \nV vs. In a plot in Fig. H. 

The features in 02,3,4(0) discussed above evolve in an 
interesting way for the models with lower noise corre- 
sponding top > 1 — see Fig. |l5|(b). The plots of 02,3,4(0) 
in this figure are produced for the model with Ao = 0.2 
and p = 3 in the same way as above for p = 1. Note that 
relative changes of the area as a function of 9 are smaller 
than for p = 1 . One reason for this is in weaker curvature 
variation for flat particles, which makes particle area less 
sensitive to the growth point position. Another reason 
is that at p > 1 the particles corners have no cusps, and 
thus particles with small areas do not appear. 



APPENDIX B: DISCUSSION OF THE 
NUMERICAL METHOD 



(i) Coarsening of the mappings which correspond to 
particles sufficiently deep in the stagnation regions. It 
was demonstrated in Ref. Jig) ] that an accurate enve- 
lope of the cluster can be obtained by using truncated 
Laurent series of F n (z). One can implement this obser- 
vation as follows. At the growth step n choose some 
1 < to < n in such a way that all particles with the num- 
bers < to are located sufficiently deep inside the stag- 
nation part of the cluster. Then one can replace the 
mapping F n = fx o ... o /„ by 

^approx) = [h Q Q /m]truncated c /m+1 o ... o fn , (Bl) 

where the mapping in parentheses which is replaced by 
truncated series is nothing but F m . One can choose to 
so that the finite series representation of the mapping 
F m (z) is accurate for z in the active growth region. By 
this trick, instead of computing a composition of n func- 
tions, one has to deal with only n — m functions at each 
growth step. Since at large n most of the particles are in 
the stagnation regions, one may have n — m <C n. 

(ii) Evaluating the particle area with lower precision. 
We used several hundred points on each particle's bound- 
ary, which produces areas accurate within 0.1%. Prac- 
tically, such a high precision may not be necessary. In- 
stead, one can predict particle areas by estimating the 
Jacobian at several points chosen within the bump ac- 
cording to some rule or randomly. 

(iii) Using an area acceptance window to discriminate 
against very small particles with areas <C a*. These par- 
ticles practically do not change the structure of the clus- 
ter branches, except near the corners between adjacent 
particles. However, due to the presence of small particles 
additional mappings appear in the composition sequence 
f n o ... o fi, which slows down the computation. 

We have not used these procedures in the simulations 
described above and neither systematically studied their 
efficiency. We felt that, at the initial stage, keeping 
growth algorithm as precise and simple as possible, even 
at the price of somewhat slowing it down, makes the re- 
sults more solid. 



Here we comment on the optimal choice of the nu- 
merical procedure. First, since the areas of new parti- 
cles are evaluated before the particles are accepted, one 
could, instead of eliminating large particles, change the 
growth algorithm so that all particles areas become equal. 
This can be achieved by adjusting the parameter X n for 
each particle until its area converges to a given value. 
Although this would somewhat slow down the speed of 
computation, an obvious gain would be in a more direct 
relation with the standard DLA growth. 

Also, one could attempt at increasing the speed and 
efficiency of the growth algorithm by introducing in it 
various improvements: 
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